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Abstract 

We develop a simple model to study the heavy quarkonium-hadron cross section in the high 
energy limit. The hadron is represented by an external electric color field (capacitor) and the 
heavy quarkonium is represented by a small color dipole. Using high energy approximations we 
compute the relevant cross sections, which are then compared with results obtained with other 
methods. 
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I. INTRODUCTION 



Quarkonium-hadron cross sections (a®h) are a necessary tool to understand the forth- 
coming data on quarkonium production, which will become available at RHIC. In the last 
six years many efforts have been devoted to this problem and real progress has been 
achieved, especially in what concerns the cross sections at low energies, close to the dissoci- 
ation threshold. In the energy region far from threshold the situation is less clear and even 
the energy dependence is still subject of debate. Extrapolation from calculations valid at 
low energies points to different directions. Results obtained with the non-relativistic quark 
model 2j indicate a rapidly falling cross section. This behavior is due to the gaussian tail of 
the quark wave functions used in the quark exchange model. This same behavior could be 
found within chiral meson Lagrangian approaches with the introduction of y/s dependent 
form factors 0|. In QCD sum rules the cross section was found to be monotonically 
increasing with energy. 

The calculations of a<$>h designed to be valid at high energies (y/s ~ 20 GeV) are quite 
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extrapolation , the model of the stochastic vacuum (MS V) 

formalism Q|. During the last years the leading order BP approach has been used very 
often. However, the recent next to leading order calculations presented in [lOj] show that, 
for the charmonium, the formalism breaks down because this system is not heavy enough. 
Most of the calculations mentioned above predict a rising cross section. In Ref. jjj, a<$>h falls 
with the energy and in Ref. Q| it stays constant. 

If the quarkonium is treated as an ordinary hadron, its cross section for interaction with 
any other ordinary hadron must increase smoothly at higher energies, in much the same way 
as the proton-proton or pion-proton cross sections. The underlying reason is the increasing 
role played by perturbative QCD dynamics and the manifestation of the partonic nature of all 
hadrons. However this partonic picture starts to be dominant only at much higher energies 
(y/s) 100 GeV). In the energy region relevant for RHIC physics non-perturbative aspects 
are still very important. In the high energy calculations mentioned above, different non- 
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Since there are still discrepancies concerning numbers (which may vary by one order of 
magnitude for different estimates) and the energy behavior, we think that it is interesting 
to calculate a$h with a non-perturbative approach, putting emphasis on the role played by 
the chromoelectric fields. In ^ | a similar treatment was adopted to study the quarkonium 
dissociation inside a QCD plasma. The color electric fields appearing in the transition matrix 
element were related to the color charge density of the medium, which, in turn, was computed 
in a specific model of the QGP. Here we start with a similar expression for the transition 
amplitude but, since we are in a purely hadronic phase, we must know the chromoelectric 
field inside nucleons and pions. There has been progress in the study of these fields, coming 
from models of the QCD vacuum Q], from lattice QCD from the Field Correlator 
Method (FCM) [17] and from Coulomb gauge QCD jig]. We hope that we can benefit from 
these advances and use the profiles of the chromoelectric fields estimated in these works 
in our problem. For this purpose, we treat the interaction between the quarkonium and 
hadron as being analogous to the interaction of a small dipole traversing a large capacitor 
and interacting with the color electric field but not with its sources. In the final part of this 
work we discuss the validity of this last assumption. Using a contact interaction between a 
heavy quark (or antiquark) and a quark (or antiquark) we compute the corresponding cross 
section and find that it is indeed much smaller than the heavy quark-external field cross 
section. The model developed here bears some resemblance to the Bhanot-Peskin picture, 
but is much simpler. Some simplifying assumptions are used to render the calculations quasi 
analytic and preserve the understanding of the basic physics. 



II. THE MODEL 



A. The interaction Hamiltonian 



The starting point is the assumption that the quarkonium (dipole) is small compared 
with the hadron (capacitor). As a consequence, the Q — Q pair will interact mostly with 
the external color field but not with the (quark) sources. Moreover, the external color field 
is considered to have only low momentum components ("soft gluons") and thus is able to 
transfer only a small amount of energy, which will be barely enough to dissociate the bound 
state. In the case of the charmonium, the typical binding energy is e ~ 0.6 GeV. Therefore, 
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in a first approximation 

e < M$ (1) 

where M$ is the mass of the bound state (M$ ~ 3 GeV) . In the case of the bottonium this 
approximation is even better. The binding energy is also small compared to the collision 
energy 

e < y/s (2) 

Inequality ((T} justifies the use of quantum mechanical perturbation theory (the Born 
approximation) and inequality (J2J) justifies the use of the eikonal approximation, which, in 
this case, implies that the hadron follows a straight line trajectory and remains essentially 
undisturbed during the interaction. In Figure 1 we present our picture of the scattering 
and our choice of coordinates, in the quarkonium rest frame: f% and r 2 are the quark and 
antiquark coordinates and E a is the chromoelectric field in the projectile, which will be a 
proton or a pion, moving with constant velocity v at impact parameter b. 

With these assumptions we can write the interaction Hamiltonian as: 

H mt = g(T?EZn + T h 2 E b 2 r 2 ) (3) 

where T a (T b ) are the generators of color group SU(3) in the fundamental (conjugate) 
representation. E% and E\ are the chromoelectric fields generated by the hadron in motion 
(capacitor) and "felt" by quark and antiquark in the bound state respectively. They have to 
be Lorentz transformed to the quarkonium rest frame, bringing to our calculation a Lorentz 
gamma factor, which is the source of the energy dependence (a/s) of our results. We shall for 
the moment neglect the magnetic component, since it does not do any work on the charges 
and thus is not effective in the energy transfer. Besides, the magnetic interaction is inversely 
proportional to the quark mass, being thus suppressed. 
We can represent this external field by: 

*(r„ t) = ,K eXp (-2LZ*!! _ exp (-^^l) (4) 

with e = 1, 2 . X, Y and Z are the hadron coordinates and 7 is the usual Lorentz factor. 
Z = vt, because the hadron moves with velocity v along the z axis. Eft, which will be 
abreviated by E, is the color electric field at the center of the projectile. The projectile 
mean square radius is related to the parameter d through: 



; r 2) = 0.86 d 



We neglect the deflection of the hadron trajectory, because we are studying reactions in the 
high energy and non-perturbative regime, i.e., with low momentum transfer. X and Y are 
related with the impact parameter b by: b 2 = X 2 + Y 2 . Notice that, by simplicity, we choose 
one preferencial direction for the field, in this case, the x-axis. 
Neglecting the CM motion, (j3J) can be rewritten as 

H m t = g{^El + ^E b 2 ){^^) (5) 

Also for the sake of simplicity, when working with (J3J), we will take x% — x 2 ~ a, where a is 
the typical separation between quark and antiquark. Initially the quark- ant iquark pair is in 
a localized region of the space. 



B. The initial state 



The initial wave function of system has spatial and color parts defined by: 

*i = f{n,r 2 )c n d n (6) 

where c n and d n , with n = 1,2, 3, are the initial color vectors for quark and antiquark 
respectively, taken in a color singlet state. We choose 

-» 9 -» 2 

T\ T 2 

f(n, r 2 ) = Ni exp[ -} exp[ -] exp(-iEit) (7) 

cr cr 

where £j (£j = M$) is the quarkonium initial energy and Ni is a normalization constant 
given by: 

The initial wave function describes the confinement of quarks and also asymptotic free- 
dom, as it allows the quarks to be independent inside the bag. It is easy to see that the 
connection between the quarkonium mean square radius and the parameter a is 



C. The final state 

Under the action of the external field the initial wave function evolves to a final state 

*f = t(r l ,r 2 )c j d k (8) 



where Cj and dk, with j, k = 1,2, 3, are the quark and antiquark final color vectors and 
t{r\, r%) the spatial part of the wave function. In the final state of this reaction we have to 
deal with the transition of a pair of an excited quark and an antiquark to a pair of mesons 
D-D (or B — B). This transition is highly non-perturbative and has to be modelled. We 
shall use here two approaches. 



We first assume that the quark and antiquark are converted into two free mesons (a M 
and a M) which are thus described by plane waves: 



where p\ and f>2 are the meson momenta and Na is a normalization constant given by: 




with V being an arbitrary normalization volume, which will be cancelled in the calculation 
of the cross section. In the above expression Ef is the final energy of the Q — Q pair. The 
energy transferred during the reaction must be sufficient to dissociate the bound state into 
a pair of mesons with open charm (DD) or beauty (BB) and therefore: 



where tum { m ji) is the mass of the meson coming out from the fragmentation of the quark 
(antiquark). With this definition of Ef we implicitly account for the conversion of quarks 
into hadrons, a process which cannot be better described in this simple model. 

The assumptions Q and (fTU|) are reasonable but they represent a case of "extreme 
freedom" : they do not take into account the energy loss from a parent quark when it is 
converted to a (less energetic) final meson. This process is described, in certain situations, 
by the fragmentation functions. Morevover, the final mesons can have any momentum 
and even though higher momenta will be naturally suppressed in the calculation, we are 
overestimating the phase space of the reaction. 



model A 



^(ri,r 2 ) = Na exp{ip\.r{) exp (ip2-r 2 ) exp(-iEft) 



(9) 




(10) 



model B 
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Given these weak points of (JHJ) and (fTUjl we shall also use a second approach for the final 
state which is more conservative. We shall assume that the energy transferred to the heavy 
quarkonium $ will transform it into an excited (but still bound) state The mass of this 
excited state will be taken to be slightly higher than the first charmonium and bottonium 
excitations and T' respectively. It is known that these excitations are very weakly bound. 
Therefore, by choosing slightly higher masses for them, which are above the D-D and B- 
B decay thresholds, we are simulating a fragmentation process to a pair of nearly at rest 
mesons. This assumption is complementary to and ()10|) since here we give to the heavy 
quarks only the "minimal freedom". The ground state wave function was chosen to be the 
Gaussian (J7j). Taking the harmonic oscillator as inspiration, we choose the wave function of 
the first excited state as a function which is odd in the x direction (this is the direction of 
the chromoelectric field) and symmetric in x\ and x 2 - 

tB(r 1 ,r 2 ) = N B — 2 exp[-—^-] exp[--^} exp(-ie f t) (11) 

where the normalization constant is: 

and a' is related to the size of the state or T'. Using the wave function (jllj) has some 
advantages. In first place, it avoids the definition of a fragmentation mechanism with the 
introduction of new parameters. In second place, as it can be seen, (JTTJ) is orthogonal to Q, 
so that the matrix element < f\H int \ty i > is zero if the Hamiltonian is a constant. Notice 
that this does not happen when we use (jHJ) and therefore the approach A might contain 
spurious contributions. The same comment is valid for the calculations made in Ref. jl^ . 
This makes the contrast between approaches A and B even more necessary. Finally, in what 
follows we shall use the Hamiltonian Q, without the approximation x\ — x 2 — a made in 
model A. 

Transition amplitudes and cross sections 

The transition amplitude for model A can be easily computed from (jHJ), ©, © and 

T fi = (V f \H int \%) = JdtJ d 3 n J d 3 f 2 ^(Fi^iM^Wifa,^) (12) 



An analogous expression holds for model B with the use of @, (@J), ©, © an d (JHJ). We 
next take the amplitude squared |T/j| 2 = T^Tfi and since color is not observed, we take the 
average of the all initial color states and the sum of all final states: 

1 ^ 1 



W^£^£££|i>| 2 (is) 

The cross section with model A is given by: 

W) Sti I W) d ° n2w L dbb ™ (14) 

The above expression is very simple and can be calculated almost analytically. Because 
of the gaussian Ansatz (J1J and © we can easily integrate (fT2|) over the coordinates and 
over the impact parameter. In the last step of (JHJ), the integration over the phase space 
had to be done numerically. In we made the additional assumption that the outgoing 
mesons are nearly at rest and we could thus simplify and perform the integration over 
Pi and f>2 analytically. Here we prefer to be more "exact" and perform the last integrations 
numerically. 

The cross section with model B is simply given by: 



a B = 2tv J dbb\T fi \ 2 (15) 
which, after the proper substitutions and integrations yields: 

aB = Y n {9Eo) (a» + a*na*a* + <P(a» + a*)]* ^ " " 2( 7 2 - 1) ) (16) 

where: 

u = e f -ei = M$/-M$ (17) 

From the above expression we can observe that the cross section rises with the energy 
(7) and saturates at a constant value. The enhancement of the chromoelectric field is tamed 
by the Lorentz contraction of the projectile. As for the size parameters, a, a' and d, the 
cross section first rises and then falls with increasing values of the parameters. The values of 
the maxima strongly depend on the model and might change for a different choice of wave 
functions. However, the physical picture is very simple. Expression (jlfi)) tells us that the 
probability of converting a quarkonium of given initial size a to a final state with size a' 
tends to zero if a' = or if a' — > 00 because the overlap between these very different states 
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and the initial state is zero. For the same reason the cross section vanishes for a = and for 
a — > oo. The parameter d is associated with the extension of the capacitor. When it goes to 
infinity the spatial dependence of the potential disappears, it becomes a constant and then 
< ^ f \H int \^i >^< V f \% >= 0. 

D. The interaction with the sources 

In the introduction it was assumed that the quarkonium is well represented by a small 
dipole, which traverses a large capacitor. However this may be a too strong assumption 
because the dipole is not always so small. For example, comparing the size of the charmonium 
with the size of pion we have tipically ^ ~ ^| ~ 0.67. Therefore it is necessay to include 
the interaction between the quark and antiquark in the quarkonium with the sources (the 
"plates" of the capacitor) which may be either a quark and an antiquark in the case of the 
pion or quark and a diquark in the case of the proton. 

In order to take these interactions into account we shall assume that the interaction 
between a quark (or diquark) in the capacitor and a charm quark (or antiquark) in the 
dipole can be divided into a short distance and a long distance part. The later was already 
included before in the interaction with the chromoelectric fields produced by the sources. 
The former will be modelled as follows. 

model C 

The short distance interaction can be approximated by the contact interaction part (the 
one with the delta function) of the one-gluon exchange potential 21\: 

H mt = Voce = £ £ %Z • - a t ■ rij* 8 ^)) (18) 

4 Kr H 3m i m J 1 

where A and a are the Gell-Mann and Pauli matrices respectively, which are responsible for 

color and spin interactions. The Coulomb term in the above expression will be neglected 

because it is of long range. The labels i = a,b and j = 1,2 refer to particles in the capacitor 

and dipole respectively. With this notation, in the interaction between particle a and 1 the 

delta function above takes the form: 

S 3 (f a - n) = 5(x a - x x ) x 5(y a - y x ) x 5{z a - z x ) (19) 

9 



where f[ = is the same as before and f* a = (x a ,ya,z a ) is the coordinate of the 

particle a in the quarkonium rest frame. In order to compute the transition amplitude we 
need to know the new wave functions, which now include both the quarkonium and the 
capacitor. They are: 



= f(n,r 2 )g(r a ,r b ) 



(20) 



and 

*/ = tc{fl, r 2 )g{r a , f h )cidjeih k (21) 

In the above expression the function / is the same as before and given by (JJJ). The function 
to represents the spatial distribution of the heavy quarks in the final state, which is assumed 
to be an excited but still bound state, very much like in model B. However, if we would 
choose to = ts, the transition amplitude < ^f\H int \^i > would vanish because the contact 
interaction does not depend on the coordinates and hence < > is the product of an 

odd by an even function of x, being thus zero. Since we are mostly interested in knowing 
the order of magnitude of this contact interaction we shall approximate the final state wave 
function by a gaussian, given by: 



tc{ri,f 2 ) = N c exp(—±) exp(—£) e 



-ie f t 



with the normalization constant given by: 



N, 



c 



1 7r ' a 
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(22) 



(23) 



The computation of the contact interaction requires the knowledge of the positions of the 
quarks in the capacitor, which is given by the function g 



9{r b , r b ) 



Np exp 



y-{x a -Xf- 








exp 


L d 2 J 





r (x b -Xf- 




\-(y b -Y) 2 ] 


exp 


\- d 2 


exp 


L d 2 J 



x exp 



-i 2 (z a -zy 



exp 



d? 



(24) 



where Z = vt, d and 7 have the same meaning as before and Np is the normalization 
constant of the projectile wave function, given by: 

,2 



Ni 



87' 

7T 3 d 6 



(25) 
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Notice that g is the same in the initial and in the final state. This assumption is consis- 
tent with the eikonal approximation introduced above and avoids the introduction of new 
parameters. 

With these ingredients we can evaluate the transition amplitude: 

T fi = (* f \H int \%) 

dt J rf 3 firf 3 r 2 J d 3 f* a d 3 f b ^}(fi,r 2 ,r^ a ,n) H in t(fi,r 2 ,r^ a ,r b ) #*(ri, r 2 , f a , f b ) (26) 

and the cross section: 

2 10 2 / 1 1 1 1x2 
a C = ~^ 7Ta s ( + + + ) 

3 4 ^m a mi m a m 2 m b mi m b mi J 

2 6/6 7 2 ° /2 n 2 i A2 

7 0_0_ /_ 2 {aP+a?) + d \ , 7) 

7 2 -l {a' 2 + a 2 ) 5 [d 2 {a' 2 + a 2 ) + 2a 2 a' 2 } P \ 4( 7 2 - 1) ' 1 ' 

where we have used f|13|) and the analogous expression for the sum and average over spins. 
Apart from a numerical factor, (fTT)|) and (|27jl have the same energy dependence. This is so 
because the same Lorentz contraction in the exponent of the Hamiltonian (jSJ) and (@J) leading 
to (jl6j) is now present in the capacitor wave function ()24|). Moreover, the same Lorentz 7 
factor, previously multiplying the E a field in Q reappears now in the normalization constant 
(|23|) . The dependence of (|2T|) on a and a' is qualitatively the same as the one found in (fTT5|) 
and has the same physical origin. Finally, the cross section above is now a monotonically 
decreasing function of d. The observed behavior with d means that, in a larger capacitor 
the quarks are spread across a larger transverse area and it becomes more difficult for them 
to find the charm quarks in the target and suffer a contact interaction. 



III. RESULTS AND DISCUSSION 



In the numerical estimates presented below, we shall adopt d = 0.8 and 0.6 fm for the 
proton and the pion respectively. We shall also take a = 0.4 and 0.2 fm for the J/ip and T 
respectively and a' = 0.8 and 0.45 fm for the ^' and T'. The bound states \l/ (m^ = 3.07 
GeV) and T (mj = 9.46 GeV) will be, in model A, dissociated into pairs of mesons D 
{m D = 1.87 GeV) and B {m B = 5.27 GeV). The excited states used in models B and C 
have masses m$/ =3.8 GeV and m$/ = 11 GeV in the case of charmonium and bottonium 
respectively. The value of the strong coupling constant and the constituent quark masses 
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are the same used in j2l|, i.e., a s = 0.64, m q = 0.3 GeV, m c = 1.2 GeV, m& = 4.74 GeV 
and the diquark mass is m<i = 0.60 GeV. 

As it is clear from (jSJ) and , we need to know the average value of the color electric field 
in the projectile gE = (h\gE\h). In a first approximation this number might be identified 
with the string tension k ~ 0.18 GeV 2 or n ~ 0.9 GeV j fm. The string tension calculated in 
[hs! ] is somewhat larger. In 17] the transverse profile of the string was studied. The strength 
of (h\gE\h) depends on the quark-antiquark (or quark-diquark) separation, being larger for 
larger systems and so far it has been calculated only for large systems. Therefore (h\gE\h) 
is another source of differences between a proton and a pion projectile. Taking an average 
of the values found in we choose (h\gE\h) = 1 GeV/fm. 

As mentioned in the introduction, our model has common aspects with the BP approach. 
Therefore we shall, in what follows, compare our results for cr$h with those obtained by 
Kharzeev in [jj: 



A 



6.5 



with A given by 



o$ fc = 2.5 I 1 - y j mb (28) 



and Ao — (M^ + e), where is the projectile mass and 

e = 2m M -M $ . (30) 

In Figure 2 we show the cross sections for the proton-charmonium dissociation obtained 
with model A (dotted lines) and model B (dashed lines) and compare them with the BP 
cross section (solid line with stars) given by (J28j) . The two upper curves are obtained with 
(h\gE\h) = 1 GeV/fm and the two lower curves with (h\gE\h) = 0.57 GeV/fm (model A) 
and (h\gE\h) = 0.53 GeV/fm (model B). With these smaller values of the chromoelectric 
field our curves come close to (|2*Hj) . Figure 3 shows the corresponding cross sections for the 
proton-bottonium dissociation. Again, the two upper curves are obtained with (h\gE\h) = 1 
GeV/fm and the two lower curves with (h\gE\h) = 0.69 GeV/fm (model A) and (h\gE\h) = 
0.49 GeV/fm (model B). As in the previous figure, reducing the value of (h\gE\h) leads to 
some agreement with ()28|) . Given the conceptual resemblance between our model and the 
BP one, it is reassuring to find a certain similarity between the results, both in magnitude 
and energy behavior, once an appropriate value of {h\gE\h) is chosen. 
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In Figure 4 we show the cross section for J/ ip dissociation by pions compared with results 
obtained the meson exchange model jjj (thin dotted line), with the quark exchange model 
P| (thin long dashed line), with short distance QCD (the BP approach) Eq. (J25j) (thick solid 
line) and QCD sum rules (thin solid line). In spite of the fact that, at such low energies 
our approach looses validity, it is, nevertheless, interesting to observe that our curve is in 
the center of the region covered by the other calculations. In Figure 5 we compare the cross 
sections p — J/ip (upper curves) and ir — J/ip (lower curves) calculated with models A 
(dotted lines) and B (dashed lines). In the high energy limit, where both cross sections are 
nearly constant, we observe that the relation between the cross sections is: 

<T P _$ ~ ScJt,-® model A (31) 

cr p _$ ~ 4.2 cr,,--* model B (32) 
which in both cases is much larger than the one expected from the additive quark model: 

o-p_$ ~ - o^-® (33) 

This is remarkable since the additive quark model relation holds for other high energy 
scattering processes like it — p and p — p. Since (h\gE\h) was kept the same for both cases, 
this unexpected relation between the cross sections must come from differences in the wave 
functions. In Figure 6 we repeat this comparison for the reactions p — T and 7r — T, finding 
(|31|) for both models. We have kept (h\gE\h) = 1 GeV/fm for both projectiles. Taking 
(p\gE\p) > (TT\gE\n) would increase the deviation from (jHSJ). 

In the high energy limit ordinary hadrons are expected to have a geometrical total cross 
section. Since the quarkonium dissociation discussed here is a more specific reaction it is 
not obvious that its cross section follows a geometrical behavior. Such a behavior was found 
in ^]: a<ph oc a s a^ where a Q is the Bohr radius of the quarkonium. In our case, as it can be 
seen from ()16|h (J27j) and from the numerical evaluation of ()14|) . we have a very non-trivial 
dependence on a. Since the initial state (containing the variable a) is the same, the difference 
between models comes from the spatial dependence of the final state. The plane waves in 
model A have no spatial scale. Therefore they are more "inclusive" and so a a should be 
closer to the quarkonium-hadron total cross section than erg. In model B the quarkonium 
ground state is converted into a resonance-like state, which wave function contains the size 
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parameter of the resonance and distorts the final geometrical behavior. Therefore model A 
is closer to a geometrical behavior than model B. 

In order to see how far we are from the geometrical behavior, we show in Figures 7 and 8 
the dependence of a a (dotted line) and <jb (dashed line) on a for charmonium (Fig. 7) and 
bottonium (Fig. 8) dissociation. The cross sections are divided by a 2 so that geometrical 
behavior translates into a horizontal line. We see that, whereas model A tends to this 
behavior, model B is far from a geometrical behavior. This indicates again that our model 
is very sensitive to the choice of the final state wave functions. 

In Figure 9 we show the cross section oc (|2*7|) for charmonium dissociation by protons 
(solid line) and by pions (dashed line). In Figure 10 we show the same quantitity for 
bottonium dissociation. We use the central values for a, a', d and a s . We can see that, 
in all processes, the cross sections are more than two orders of magnitude smaller than the 
corresponding cross sections computed with model A or model B. No possible change in 
parameters could make these cross sections comparable. Another feature of these curves is 
that the cross sections for J/ip dissociation by pions are larger than those for protons by a 
factor close to 4. This might be guessed looking at ()27jl . The pion is light quark- ant iquark 
system and the proton is light quark-diquark dipole. The diquark is twice heavier than a 
constituent quark. Whereas for the pion m& = m a , for the proton we have m& = 2m a . 

Before concluding we would like to make a remark concerning medium effects on the 
cross sections calculated above. We are primarily studying reactions which happen before 
thermalization (in nucleus-nucleus collisions) or with no thermalization at all (in proton- 
nucleus collisions). The formation time of the heavy quark pair is of the order of 0.2 fm. The 
thermalization time of hadronic matter formed in heavy ion collisions is a model dependent 
quantity. Early estimates pointed to 1 fm. Recent estimates point to 0.6 fm. Even taking 
seriously this last number, it is fair to say that heavy quark pair production (and collision 
with a hadron at high energies) precedes the formation of an equilibrated medium. After 
thermalization, the energy is completely redistributed and collisions occur at energies of the 
order of the temperature (< 1 GeV). In this regime we do not expect our approach to be 
valid. The effects of a thermal medium on the heavy quarkonium are known |24j: the string 
tension becomes weaker, the quarkonium size increases and its mass decreases. These effects 
are, all of them, very small except if we get close to the deconfinement transition temperature. 
In view of these considerations we have neglected medium effects in our calculations. 
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IV. CONCLUSIONS 



We have developed a simple model for the non-perturbative quarkonium-hadron interac- 
tion. At the present stage of the field, this sort of model is still useful to organize the ideas. 
We tried to make simple and yet realistic choices for the interaction Hamiltonian and for 
the wave functions. In particular we have treated the final state in two very different and 
complementary ways. Simple models are not appropriate to provide very precise results but 
they can help in determining the order of magnitude of the cross sections and their behavior 
with the reaction energy. Having said that, we can summarize our conclusions as follows. 

i) The charmonium-hadron cross section is of a few milibarns. The bottonium-hadron cross 
section is about four times smaller. This is in agreement with most of the previous calcula- 
tions. 

ii) All cross sections grow with the reaction energy and reach a plateau in the high energy 
limit. This is in agreement with the BP approach. 

iii) In this limit they do not obey the simple relations derived from the additive quark model. 

iv) Also in this limit our cross sections deviate significantly from the geometrical behavior 
(a oc a 2 ). 

v) The contact interactions between the heavy quarks and the light quarks in the light 
hadrons is negligible compared to the long-distance quark - E a field interaction. This is sur- 
prising since sometimes the dipole and the capacitor have similar sizes. This finding gives a 
posteriori support to our model and also to the BP approach. 

Conclusion i) may be relevant for RHIC and LHC physics. Conclusions ii, iii and iv 
suggest that the heavy quarkonium has interaction properties which are very different from 
light hadrons. This has been conjectured before. In particular, in 22] this difference was 



attributed to the fact that in heavy quarkonia the energy is mainly stored in the masses 
whereas in light hadrons the energy (mass of the hadron) comes mostly from the gluonic 
fields. 
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FIG. 1: Quarkonium-hadron interaction in the quarkonium rest frame. The hadron is a capacitor 
moving to the right. The E a field is along the x direction. The thick black dots are the capacitor 
"plates": quark-antiquark for the pion and quark-diquark for the proton. 
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FIG. 2: J/ip — p cross section with model A (dotted lines), model B (dashed lines) and with the 
Bhanot-Peskin approach (lines with stars). Upper curves: stronger E a field. Lower curves: weaker 
E a field. 
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FIG. 3: Same as Fig. 2 for the T + p cross section. 
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FIG. 4: Pion-charmonium cross section as a function of y/s with several models. 
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FIG. 5: Charmonium-hadron cros section with model A (dotted lines) and model B (dashed lines). 
Upper curves: ajj^ — p. Lower curves Oj/^ — n. 
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FIG. 6: Same as Fig. 5 for bottonium-hadron cross sections. 
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FIG. 7: Charmonium-hadron cross sections as a function of the charmonium size parameter for 



model A (dotted line) and for model B (dashed line). 
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FIG. 8: Same as Fig. 7 for bottonium-hadron cross sections. 
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FIG. 9: Charmonium-proton (dashed line) and charmonium-pion (solid line) cross sections calcu- 
lated with model C (contact interactions). 
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FIG. 10: Same as Fig. 9 for bottonium-hadron cross sections. 
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